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ABSTRACT 

Strong magnetic fields modify particle motion in the curved space-time of 
spinning black holes and change the stability conditions of circular orbits. We 
study conditions for magnetocentrifugal jet launching from accretion disks around 
black holes, whereby large scale black hole lines anchored in the disk may fling 
tenuous coronal gas outward. For a Schwarzschild black hole, magnetocentrifu- 
gal launching requires that the poloidal component of magnetic fields makes an 
angle less than 60° to the outward direction at the disk surface, similar to the 
Newtonian case. For prograde rotating disks around Kerr black holes, this angle 
increases and becomes 90° for footpoints anchored to the disk near the horizon of 
a critically spinning a = M black hole. Thus, a disk around a critically spinning 
black hole may centrifugally launch a jet even along the rotation axis. 



Introduction 



The production of relativistic jets in (micro) quasars, Active Galactic Nuclei and Gamma 
Ray Bursts involves the conversion of rotational energy of the accreting matter and/or of 



the spinning black hole into plasma bulk motion (Blandford & Znajek 1977 Blandford & 



Payne 1982 Punsly & Coroniti 1990). This energy conversion is mediated by the magnetic 



field, which can both accelerate plasma by magnetic gradient forces and collimate it by 



hoop stresses, (see Meier et al. 2001, for review). It is not clear whether jet is powered by 
the accretion disk or by the spinning black hole (more precisely, whether the footpoints of 
magnetic field lines that eventually form the jet should pass throughout the disk, ergosphere 
or horizon of the black hole). Numerical simulations confirm, at least partially, that all 



these mechanisms can be operational (Meier et al. 2001 Koide 2004 Komissarov 2005 



Tchekhovskoy et al. 2008) 



There are three related issues in the production of jets: launching, acceleration and 
collimation. In this paper we address the conditions for jet launching from a thin accretion 
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disk rotating in the equatorial plane of a Kerr black hole. We assume that the disk is 
permeated by a large scale magnetic field, so that magnetic field lines are anchored in the 
disk and are dragged both by plasma rotation and the rotation of space time around the 
black hole. This assumes that the disk conductivity is very high, and that inside the disk the 
magnetic field is subdominant and does not affect plasma motion in the disk. In addition, 
we assume that the disk is surrounded by a tenuous coronal plasma and the magnetic field 
is sufficiently strong to affect particle motion in the corona. The dynamical effects of a large 
scale magnetic field on single particle motion are often approximated as a guiding wire with 
a particle playing the role of a bead. This simple approximation restricts particle motion 
across the field, neglecting various cross-field drifts. 

We study the dynamical behavior of coronal particles, located close to the disk, which are 
locally restricted to move exclusively along a given magnetic field line. In such a case, under 
certain conditions discussed below, a coronal particle executing a circular motion around a 
black hole is in a state of unstable equilibrium, even though without a magnetic field it may 
be in a stable equilibrium. The condition of unstable equilibrium, we assume, constitutes a 
condition for jet launching from an accretion disk. This may be the most relevant condition 



for the generation of jets (Blandford & Payne 1982; Ogilvie & Livio 2001 Spruit 2008) 



2. Bead on a wire rotating in Kerr space 

2.1. Rotating Kerr metric 

Assume that a footpoint of a magnetic field line is anchored in accreting matter and 
is rotating with angular velocity uo at r = r . Next, neglect the toroidal component of the 
magnetic field (see Appendix [AJ and assume that a field line makes an angle 0i with the 
radial direction, Fig. [T] 

As a mathematical problem, consider the motion of a bead on a wire sticking out of the 
equatorial plane of a Kerr black hole at an angle $i with the radial direction and confined to 



move in a given meridional plane. In Boyer-Lindquist coordinates (e.g. Misner et al. 1973), 
make a coordinate transformation to the frame rotating with the foot-point, 0^0' — uot. 
Restricting particle motion along the wire, tan# = r/((r — ro)tan6 1 ;), dcf) = (see Fig. [T]) 
implies 

d» = oot« I? J°f 2 (1) 

(r — r ) 2 + r 2 cot 2 &i 

In the coordinate frame rotating with the wire the non-vanishing components of the metric 
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Fig. 1. — Geometry of the model. A magnetic field line is anchored in the disk at the 
equatorial plane of a Kerr black hole at radius rn and makes an angle Q\ (in Boyer-Lindquist 
coordinates) with the radial direction. 
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where a is the angular momentum per unit black hole mass and M is the mass of the black 
hole. (We set G = c = 1, use (—1, 1, 1, 1) sign convention and assume that the mass of a 
test particle is unity.) We assume that uo > 0, while the parameter a can be in the range 
— M < a < M (a > corresponds to prograde rotation of the wire). Below we refer to the 
metric ^ SIS db rotating Kerr metric. 

Condition g 00 = defines two light cylinders. Let us consider the location of the light 
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cylinders for a particular case of radial magnetic field Q[ = 0, 9 = n/2. In this case 



goo = -1 ~ (a 2 + r 2 )u 2 + 



,j. 2 



2M(1 - acj) 2 



9rr = -£- (3) 

For a; = 0, the inner light cylinder coincides with the ergosphere, r = 2M. For small 
u; << 1/M, the outer light cylinder is approximately at r ~ l/o;. The location of the light 
cylinders are given by 

3^(1 -aufMu 








sin 


-arccos 




v 3 \ 




(4) 

Since the determinant of the metric tensor is smaller than zero beyond the light cylinders, 
the approximation of a rigidly rotating wire is inapplicable in those regions. 

Condition d r goo = defines circular orbits 

rlu 2 = M(l - aujf (5) 

This is Kepler's law in a Kerr metric. 

Transformation to the rotating Kerr metric is physical only for angular velocities smaller 
than u p h, angular velocity of a photon orbit, defined by the conditions of circular rotation 



with the speed of light g o — 0, d r g o = 0. This gives — 4a 2 M + r(— 3M + r) 2 = (Bardeen 



et al. 1972). For a given a and r, the transformation to the rotating Kerr metric becomes 



meaningless for to higher than the angular velocity of a photon circular orbit, 

Wpb =p|a| + 6Mcos(|arccos(=F|a|/M) 

The upper sign corresponds to prograde rotation. Particular values are a = 0, r p h = 
3M, uj ph = 1/(3 V / 3M) for Schwarzschild black hole, a — M, r ph = M, u ph = 1/(2M) for pro- 
grade and a = —M, r ph = 4M, = 1/(7M) for retrograde photon orbits. For 1/(7M) < 
uj < 1/(2M) this requires sufficiently high a, satisfying uj < 1/ [a + 2y/2M (l + cos (|arccos(— ^)))) 
(for a = this requires uj < l/(3v3)). For to — > 1/(2M) both light cylinders merge on the 
black hole horizon at a = M = r. For higher u, the transformation to the rotating frame 
becomes meaningless everywhere (Fig. [2]). 



3/2 
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Fig. 2. — The locations of light cylinders for different uj as function of a. The light cylinders 
always lie outside the horizon (dotted curve). For a given uj, the inner and outer light 
cylinders merge at the photon circular orbit (dashed line). The physically meaningful region 
lies in between the light cylinders, to the right of the light cylinders curve. For uj — > 0, the 
inner light cylinder coincides with the ergosphere r = 2M. As uj increases, the outer light 
cylinder moves to smaller r; the radial location of the inner light cylinder is a complicated 
function of uj. For uj > 1/(7M) there is a region for sufficiently small a, for which the 
transformation to the rotating frame is unphysical. For high uj — > 1/(2M) and high a > 
M/y/2, the inner light cylinder moves inside the ergosphere of a Kerr black hole. For uj > 
1/(2M), the transformation to rotating frame is unphysical everywhere. 

2.2. Motion of particle in rotating Kerr metric. 

Since metric O) is time independent, the Hamilton - Jacobi equation, 



g^d^S = g 00 (d t S) 2 + g rr {d r Sf = -1 



(7) 
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has a solution in the form S = —Ent + S r (r) with 



S r = J dr^-(E 2 + g 00 ) 



Qrr 

9oo 



Differentiating with respect to E Q gives the equation of motion 



(9 t r) 2 = -(l + W^o)— (9) 

9rr 

Transforming to proper time dt/dr = E /(—g 00 ) gives 

(d T r) 2 = -^±^ = V (10) 

QOOQrr 

This defines the effective potential V. The extremum of the potential d r V = defines the 
angular velocity of a circular orbit. Direct calculations give, naturally, condition (J5]). 



2.2.1. Reversal of centrifugal acceleration at photon circular orbit 

The above relation simply in case of a radial magnetic field, 9i = 0. In this limit 
A 



V = — 



1 - 



El 



d 2 r 



dr 2 



1 



2A1 



d r V 



(1 — aco) 2 — u 2 (r 2 + a 2 ) 



(for circular orbits) 



Mr - a 2 fa + (r 2 + a 2 )u\ a(M -r)+u (r 2 (r - 3M) + a 2 (r + 



+ 



Eo=V-9oo 



For a Schwarzschild black hole, a = 0, Eq. (11 ) shows the reversal of centrifugal acceleration 

(12) 



at the photon circular orbit r = 3M (Abramowicz & Prasanna 1990): 

d 2 r 



Or 



4 + (r-3My 



For a Kerr black hole, no clear separation can be made between the effects of the wire 
rotation and the rotation of space-time, so that the notion of a centrifugal force becomes 



ill-defined (see 


Abramowicz & Prasanna 1990 


de Felice 


1991 


Abramowicz & Lasota 


1997 


Iyer & Prasanna 


1993 


for discussion). 



3. Stability of circular orbits in magnetic field 

Circular orbits are defined by d T r = (V = 0, E = ^—goo) and d 2 r = (V^ = 0, 
condition (JsT) ) . The stability of an orbit depends on the second derivative of the potential 
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V, Eq. ph. 



-V" 



A (3a 2 - 4a v / iWr + r 2 (l - 3 cot 2 0,] 



-M, 



(13) 



Eo = V ^ M ,r=r (r 4 (r - 3M) + 2av / Mr 7 /2)( A + r 2 cot 2 

where we replaced r by r. Stability requires the epicyclic frequency to be real, k 2 > 0. We 
identify parameters for which n 2 < as a necessary condition for jet launching. 



The sign of n 2 in Eq. (13) determines whether a particle confined to move along the 
magnetic field anchored at radius r to the disk is in a state of stable equilibrium (k 2 > 0). In 



Eq. ( 13 ), the denominator changes sign at the photon circular orbit r 4 (r— 



0, Eq. (6|. The physically realizable case corresponds to radii larger than the radius of the 



photon circular orbit. The numerator in Eq. (13) changes sign at 



cot 61 



4a 



M a 

2t~3/2 f 



(14) 



When the magnetic field lines make an angle with the radial direction less than the one given 



by Eq. (14), a particle tied to the field line in the equatorial plane is in a state of unstable 



equilibrium and will be flung away from the disk. 



An analysis of Eq. (14) reveals, first, that the magnetic field cannot stabilize particle 



motion for an arbitrary (especially radial) direction of the magnetic field (angle 81 is never 
zero). In the Schwarzschild metric, cot Q\ = a/1/3, 61 = n/3, which is exactly the same result 
as in Newtonian mechanics (Blandford & Payne 1982[ ). For prograde rotation, a > 0, the 
instability region is larger, Q\ > tt/3, so that the flow is more unstable in the Kerr metric 
if compared with the Newtonian and Schwarzschild cases. Or, equivalently, a flow my be 
launched over a wider range of angles. 

Physically meaningful cases correspond to footpoint motion along stable circular orbits 
in a Kerr space-time. For prograde rotation around a critically spinning Kerr black hole, 
the outer and inner light cylinders, photon circular orbit and horizon all coincide at r = M 



in Boyer-Lindquist coordinates; in proper radial coordinates they remain distinct (Bardeen 
et al. 1972). According to Eq. (14), in the limit a — > M, the angle 9i — > ir/2, so that 



the motion along all field lines, including those aligned with the axis of rotation, becomes 
unstable. Thus, in the case of a critically spinning Kerr black hole, a jet may be launched 
by magnetocentrifugal mechanism along the rotation axis of a black hole. 

For any a 7^ M, the condition 61 = tt/2 is satisfied only inside the horizon, so that for 
a M there is always a set of magnetic field lines around the direction of hole's rotation 
where a particle is in a state of stable equilibrium. A jet won't be launched along those field 
lines. 
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Fig. 3. — The maximum unstable angle Q\ of the magnetic field inclination to the plane of 
the disk evaluated at the innermost stable orbit . For angles less that a particle on a 
circular orbit is in an unstable equilibrium and will be flung away by centrifugal forces. For 
a Schwarzschild black hole, a = 0, $i = 60°, similar to the Newtonian case. For a — > M, a 
particle can be accelerated along the rotation axis. For counter rotating disk around critically 
spinning black hole the critical angle decreases to cotOi = Ay/2/9, 9i = 57.85° 



, the angle Q\ is, generally, close to 60°, 
except in a narrow region of parameter a ~ M for prograde rotation, where 9i approaches 
90°, see Fig. [3j For a retrograde disk rotating around a critically spinning black hole, 
the maximum launching angle decreases to cot #2 = 4y / 2/9, Q\ = 57.8° at the last retrograde 
stable orbit corresponding to r = 9M. In Fig. [4] we plot the values of the maximum unstable 
angle Q\ of the magnetic field inclination to the plane of the disk as function of radius various 
values of black hole spin a. 



At the innermost stable orbit (Bardeen et al.|1972 
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Fig. 4. — Values of the maximum unstable angle Q\ of the magnetic field inclination to the 
plane of the disk as function of radius for various values of black hole spin a. Dashed line 
is the values of Bi at the innermost stable orbit, see Fig. |3j Dotted line corresponds to 
Schwarzschild black hole, a = 0, Q\ = 60° 



4. Force-free magnetospheres 



Conditions for stability of particle motion at the base of the jet derived in ^3] should be 
applied to the magnetic structures of black hole magnetospheres to determine if a particular 
magnetic field solution is consistent with the jet launching conditions at the disk surface. 
Analytical models of black hole magnetospheres are usually limited to the force-free approx- 
imation, when matter inertia and pressure forces are neglected. The equation governing 



force-free magnetic configurations in the Kerr metric was derived by Blandford fc Znajek 
(1977) (see also Thorne et al. ( 1986[ )). No analytical solution for the magnetic field structure 
of fast spinning black holes, a ~ M, is known. (The vacuum Wald solution ( |Waid|p74l ) has 
large parallel electric fields; this limits its astrophysical applicability. In addition, the Wald 
field is orthogonal to the equatorial plane.) 
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Analytical solutions of force-free structures for Schwarzschild and slowly spinning Kerr 
black hole magnetospheres (a — ► 0) are generally limited to the case of zero poloidal current]^] 
In this case (a = 0) the governing Grad-Shafranov equation ( Shafranov|1966 ) for the poloidal 



flux function P becomes (Blandford & Znajek 1977 Thorne et al. 1986) 



r 2 d r ((1 - 2M/r)d r P) + (1 - /i 2 )^P = 



(15) 



where \i = cos 6. Magnetic fields measured by a locally stationary observer are then given 



by (Thorne et al. 1986) 



Bj 



d e P 
r 2 sin 9 



1 - 2M/r 



d r P 
r sin 8 



(16) 



(Note, that a a form of the flux function r{9) as measured by a set of locally stationary 
observer is given by dr / d6 
the flat space). 



BqJ B f 



'l — 2M /r(rd r P/dgP), not just P = constant as in 



There is a number of known solutions of Eq. ( 15 ) . Those connecting to the accretion disk 
(as opposed exclusively to the black hole) are (i) a Schwarzschild black hole in the constant 



magnetic field P oc r 2 sin 2 9; (ii) the non-separable parabolic solution of Blandford & Znajek 
(1977) P oc (r — 2M)(1— fi)— 2M(l+/i) ln(l+/x).; (hi) the separable Schwarzschild paraboloid 



P oc (r + 2Mln(r — 2M)) (1 — /x) (and corresponding higher order multipoles Ghosh 2000). 



These are all linear solutions, so that any combination is a possible solution as well. 
The angle 8b that magnetic field lines make with the equatorial plane is 



tan 6 h 



Be 

B; 



y/1 - 2M/r 



rd r P 



d ff P 



(17) 



=7r/2 



The constant magnetic field solution (i) has a magnetic field orthogonal to the disk at all 



points. For the parabolic solution of Blandford & Znajek (1977) (ii), the angle that the 



magnetic field lines make with the equatorial plane decreases as the footpoint approaches 
the horizon, tan#b = yl — 2M/r; for large radii, r ^> M, field lines stick out at 45° 



to the plane of the disk. Thus, for the parabolic solution of Blandford & Znajek (1977), 



the particles confined to move along field lines are always in an unstable equilibrium, so 
that the jet can be launched from any point on the disk. For the separable Schwarzschild 
paraboloid solution (iii), the angle that magnetic field lines make with the equatorial plane 



-'^One may also find spheromak-type solutions for linear dependence of the poloidal current on the poloidal 
magnetic flux, but they do not seem to correspond to any physically interesting case. 
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is tan 6^ = U/i - 2M/r(l + (2M/r) ln(r - 2M))j . It reaches 61 = vr/2 somewhere close 
to the horizon (at the root of 1 + (2M/r)ln(r — 2M)) = 0) and then again at the horizon 
r = 2M, see Fig. (l5j). Thus, the force-free magnetic field structures of Schwarzschild black 
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Fig. 5. — Values of the angle 6b between magnetic field and the equatorial plane for force- 
free magnetospheres of Schwarzschild black hole: Blandford & Znajek (1977) and separable 
Schwarzschild paraboloid (Ghosh 2000). These solutions, and a solution for a black hole in 
constant magnetic field, are linear solutions, so that any combination is a possible solution 
as well. A given solution can centrifugally launch a jet if the angle 6\, is smaller than the 
critical angle at the innermost stable orbit. 



holes allow for a variety of field configurations with different inclinations of magnetic field 
lines at the disk surface. We may only speculate that Kerr holes have a magnetic field 
structure with various magnetic field inclinations at the disk surface as well. 
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Conclusion 



In this paper we discussed the conditions for the centrifugal launching of jets from a thin 
accretion disk around spinning black holes. Typically, it is believed that "for acceleration 
by the centrifugal mechanism to be effective, the field lines have to be inclined outward" 



(Spruit 2008). This leads to a collimation problem: a flow should first accelerate by mag- 



netocentrifugal forces outwards in cylindrical radius and later, beyond the Alfven surface, 
it needs to be collimated along the axial direction by hoop stresses of bent back magnetic 
field lines ( Blandford fc Payne|1982 ). For relativistic motion, such collimation becomes very 



ineffective (Eichler 1993). 



We found that Schwarzschild black holes conform to these expectations: conditions for 
jet launching (defined as an instability of a particle on a circular orbit confined to move along 
a magnetic field line) remain the same as for Newtonian disks: a field should be inclined to 
the disk surface at an angle 8i < 60° (when viewed in Boyer-Lindquist coordinates) at any 
radius in order to launch a jet. On the other hand, in the case of a critically spinning Kerr 
black hole, the centrifugal acceleration along magnetic field lines anchored in the disk may 
occur virtually along the rotation axis of the black hole. Qualitatively, this removes the need 
for a flow to expand sideways to get accelerated. 

The present result has implications for the transport of angular momentum by large 
scale magnetic fields. Since the particle orbits in the inner regions of accretion disks around 
spinning black holes are unstable to large scale magnetic fields with a wide range of angles 
at the disk surface, this make the transport angular momentum away from the disk more 
efficient than in the case of Schwarzschild black holes. 

The centrifugal acceleration mechanism (a fling) isn't the only way to magnetically 
launch a jet. Alternatively, a jet may be launched by stresses of rotating twisted magnetic 
field(a spring), (see, e.g. Tchekhovskoy et al. 2008). In a "spring"-type launching the disk 



or the central black hole act as a Faraday disk, which launches an outflow with the terminal 
speed equal, approximately, to the Alfven velocity in the corona. If the coronal plasma is 
strongly magnetized, the terminal velocity may be relativistic. 

It is tempting to relate the possibility of centrifugal acceleration along the rotation axis 
in critically spinning black holes to the results of numerical simulations of jet launching. It 
is well established that the efficiency of jet production is a strong function of the black hole 



spin (Meier et al. 2001); for a rapidly spinning black hole the efficiency of converting the 



accreting energy into jet power reaches the radiation efficiency of thin disks (De Villiers et al. 



2003, 2005). We speculate that this is because jet launching for critically spinning black hole 



is more efficient due to effect discussed in this paper. 
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It is expected that for relativistic flow with rotation rates nearing the speed of light, the 
toroidal component of the magnetic field is not negligible. In this Appendix we show that 
inclusion of the toroidal component does not change the stability conditions described above. 
First note, that in the case of Keplerian (non-realtivistic) rotation, the azimuthal component 
of the wire direction does not change the stability condition, as can be trivially shown (both 
centrifugal force and the force of gravitational attraction do not have 0-component, and thus 
their projections on any oblique direction become equal for the same parameters, as in the 
case of no toroidal component of the wire). 

In Kerr metric, we can repeat the previous derivation allowing for finite displacement 
in the azimuthal direction, along a direction given by angle <pi with respect to the radial 
direction. For a given redial displacement dr the change in the azimuthal angle is 



The non-zero components of the metric tensor are goo-, g rr , g rt (not given here explicitly). It 
is important that in the rotating Kerr coordinates, inclined with respect to radial direction, 
the goo component of the metric tensor remains the same, Eq. (|2]). 

Solution of the Hamilton- Jacobi equations are now 



A. Effects of toroidal component of magnetic field 



d(j) = cot 6i 




rodr 



(Al) 




(A2) 



Equation of motion is 




(A3) 



This preprint was prepared with the A AS IATgX macros v5.2. 
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Transformation to proper time 



dr 
dt ' 

gives the effective potential 



rr 9 



~(E 2 o + 



rt 



— Eq a/ goog rr — g 



(A4) 



dr 
dr 



El + goo 



9w9r 



9rt 



V 



(A5) 



Conditions of circular orbit give V = (E = —^—goo), K' ^ 9oo,r — 0. Thus, similar to 
the case of radial wire, the epicyclic frequency of an oblique wire is proportional to goo,rr- It 
changes sign at the same vertical angle 9i as in the case <pi = 0. Eq. (14). 



B. Pseudo-potentials 

General relativistic effects are often modeled using pseudo-Newtonian potentials. For 



Schwarzschild black holes such a potential was first proposed by Paczynsky & Wiita (1980), 
&pw — M/(r — 2M). For the Paczyhski- Wiita potential and centrifugal potential (l/2)u 2 r 2 , 
the angular velocity corresponding to circular orbits is lu = ^M/(r (r - 2M) 2 ) and the 
critical angle 

t a n 2 9 PW = 3 + ^^- (Bl) 
r — 2M 

This value is close to the exact value in Schwarzschild geometry. Even at the last stable 
orbit (r = 6M, tan 2 6pw = 4) it differs from 60° by only 3.43° degrees. Thus, in the case 
of Schwarzschild black hole the Paczyhski- Wiita potential reproduces reasonably well the 
launching conditions at the disk surface. 

In principle, general relativistic effects should also be taken into account in calculation of 



centrifugal acceleration. Using relation (12) for the centrifugal force, the pseudo-centrifugal 
potential becomes (l/2)u; 2 (r — 3M) 2 ; the angular velocity corresponding to circular orbits is 
then oo = \J M / ((tq — 3M)(ro — 2M) 2 and the critical angle 

tan 2 d PW = 3 + — + — B2 

r — 2M r — 3.1/ 



At the last stable orbit (assuming it is at r = 6M) tan 2 6pw = 5, again reasonably close to 
the exact value tan 2 8 = 3. 



